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Abstract 

A new model for the intermediate compound of the hydrogen evolution 
reaction (HER) is proposed, for the electrochemical reduction of hydrogen in 
the presence of bisulfate on platinum(lll). The formation of this compound, 
a regular 2 dimensional honeycomb ice lattice, occurs by a first order phase 
transitions that involves the reorientation of water molecules. The model is 
analyzed using new and simple effective cluster approach which highlights 
the relevant transitions in system. This method is based on the cluster vari- 
ation method used successfully in our previous work on the UPD of Cu onto 
Au(lll), and permits us to explore a large region of parameter space, an 
essential feature to study this complex system. The theory makes full use of 
the properties of the diffuse layer: The water molecule is reoriented as the 
potential is changed. For positive potentials it forms linear chains which are 
responsible for the a/3 x structure of the sulfate observed by STM. At neg- 
ative potentials water turns so that its dipole points towards the Pt. Then it 
will form a regular honeycomb network of hydrogen bonded molecules, with 
the sulfate at the center of the hexagons. Then the bisulfate is desorbed, 
leaving the honeycomb HER structure behind. Our model thus provides an 
explanation of the well known fact that only 2/3 of the Pt atoms participate 
in the electroreduction of the hydrogen. The theory implies geometrical con- 
straints to the water potential: It should 

• be of tetrahedral coordination. 

• be analytical to be able to include the double layer effects, the ions and 
external fields, in the calculations. 
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• The angular part of the potential must be relatively soft, since the angle 
0-H-O of the hydrogen bond will be bent to keep the water near the metal 
surface. 

The analytical tetrahedral Yukagua model of water satisfies all of these re- 
quirements. 



1 Introduction 

The underpotential deposition of hydrogen on platinum(lll) in aqueous sul- 
furic acid was studied by Clavilier et al. using his pioneering techniques of 
preparing mono crystalline electrode surfaces []I|, 0]. 

In Fig.l we show a recent voltammogram for this system obtained by 
V. Climent on a rather perfect single crystal. For our discussion we will 
divide this voltammogram in three regions: 

1. The H-UPD region: A new structure of the HER intermediate is pro- 
posed, which explains the observed 2/3 of a monolayer yield for this 
surface. 

2. the CP butterfly region: Here there are two distinct features, a sharp 
spike, which we identify with the formation of a a/3 x a/3 lattice of 
water-like molecules, with its dipole pointing up. 

3. the positive region: In this region the a/3 x a/7 structure is formed by 
'polywater' chains of water molecules with its dipole pointing down 

We note that these steps are similar to the ones proposed by us for the 
UPD of Cu/Au(lll)[|, I, |: In reverse order: 

• Region 3: The formation of stacked, elongated hexagonal structures 
corresponding |0, ||, P] to water-bisulfate coadsorption is explained by 
the formation of a hydrogen bonded water chain[0, |ll|. There is no 
phase transition in this case. 

• Region 2: The spike here is due to the co-adsorption of bisulfate and wa- 
ter, forming in fact the same honeycomb structure as in the Cu/Au(lll) 
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Figure 1: The voltammogram for the H2/ Pt{lll) / SO4H2 (cour- 
tesy of V. Chment) system is divided in three regions: l)The 
UPD region, 2) the CP butterfly region and 3) the positive re- 
gion 
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case. Now the sulfate simply leaves an empty site, while in the copper 
case it was replaced by a copper via a phase transition. 

• Region 1: In this region the transition is different because as the sulfate 
is desorbed, its site was filled by another atom of copper. In our present 
case, the hydrogen bond of the hexagons are saturated, and the sulfate 
is simply desorbed. 

In section 2 we outline the basic formalism and our method of computing 
the fugacities of our active compounds. In section 3 we discuss the positive 
Region 3 of the voltammogram. In section 4 we present our model for the 
CP butterfly region, and finally in section 5 we discuss our model HER in- 
termediate. In the final section we compare our resulting voltammogram to 
experiment. 



2 BASIC FORMALISM 



We use a microsopic model, based on the sticky site model (SSM) |T2|, [T3|, [14 
in which the exact hamiltonian of the entire interface is projected onto an 
equivalent 'electrode' surface lattice: This means that all solvation and dou- 
ble layer effects are included. For example, for dilute solutions such as the 
one of our present example the Gouy-Chapman approximation is quite ac- 
curate, as has been shown in previous work 0. 

One of the problems in surface electrochemistry is the complexity of even 
seemingly simple systems, like the one of this paper. Even with modern day 
resources it is very difficult, if not outright impossible, to explore the pa- 
rameter space of this problem by computer simulation. For this reason we 
have used the 'Cluster Variation Method (CVM)|5|, |16|, |T|, |2| to 



calculate the adsorption isotherms and generated a voltammogram for the 
deposition of H on Pt(lll) in the presence of 0.1 M aqueous sulfuric acid. 
In our our previous work on the UPD of Cu on Au(lll)[^ |, ^, ^ we used 



the " cluster approximation" |^ . This method successfully predicted the lo- 
cal structure and transformations of the voltage dependent phases that were 
later observed. For the UPD of hydrogen and the bisulfate anion we use first 
a simplified version of this method, the PCVM ( PADE-CVM) approxima- 
tion which further simplifies the calculations. This work will be described in 
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a separate paper PT| In our present work we will use a still simpler theory 



that we call the assymptotic cluster model (ACM), where the behavior of the 
model is described by a few assymptotically large diagrams which determine 
its transitions. The interaction parameters are obtained from theory or sim- 
ulation. 



The adsorption of species i on the platinum surface primarily depend on 
the activities of these adsorbates and the lateral interactions existing between 
them|ll2|, [l3| , |T^. The activity of a species i is given by the product between a 



sticky parameter Xi{i/j) and the species contact density at the surface p°(0, 



as 



^, = \(v^)p°(o,V') . (1) 

The sticky parameter Xi{ip) is a phenomenological parameter characterizing 
the quantum mechanical interaction between the adsorbate and surface, and 



will in general, be potential dependent. We write 20 



where e is the elementary charge, Q is the partial charge of the adsorbate 
at the surface, and binding constant characterizing the overlap of 

electron orbitals between the adsorbate and surface atoms. The potential 
dependence of the contact density is given to a very good approximation by 
the Gouy-Chapman expression 

p°(0,V^)=p°(0,0)e-'^^^^('^-^-) , (3) 



where p°(0, ippzc) is the density at a smooth (i.e. without sticky sites) surface [ll4 
with a potential of zero charge, and z/j is the ionic charge of the adsorbate. 
It has been verified in previous work that for our relatively dilute concen- 
trations of sulfate ions this is a quite good approximation. 

The potential dependence of the activity Zi is then found to be 

= A,(0)p°(0,0)e-^^>'=(^-'^°) , (4) 
where 7i = z^j + Hid is an effective electrovalence of species i. 
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We obtain analytical expressions representing clusters of interacting ad- 
sorbates which arc used to calculate the coverages of hydrogen and the bisul- 
fate anion as a function of potential. The crucial parameters in our present 
theory are the water orientation and the bisulfate adsorption. 

2.1 THE ADSORPTION AND ORIENTATION OF 
WATER AND (BI) SULFATE 

We discuss first the potential dependent adsorption fugacities of bisulfate 
and water. In our model water is adsorbed on the atop position when its 
dipole is pointing down and in the hollow site when the dipole is up. 



2.1.1 BISULFATE 

The inner layer equivalent fugacity zs for the adsorption of the bisulfate is 

^5 = A°p°(0,0)e^^-/^^(^-'^r) (5) 

where (3=1/ kT is the Boltzmann thermal factor, the electrosorption valency 
of the bisulfate is ^5 = — 1, and ipg'^ is the electrosorption reference potential, 
that depends on the nature of the substrate. The sticking coefficient can be 
interpreted as XsW = exp[/3/is], with /is as the free energy change that 
occurs when a bisulfate ion binds to the metal surface. ^^(0, ip) is the inner 
layer local density of bisulfate for a local potential ip, which is estimated from 
the Gouy- Chapman formula. 

zs = A°p° (0, 0)e^<^-^^^) = zy^^^^-^sn (6) 

where we take 

^l,P'^ = 0.8V{RHE) (7) 
and the bisulfate fugacity zs is {zg = 1, T = 298. 16/^) 

^^^g38.922[^-<-]^ (8) 

The electrovalence is taken as = — 1. 
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2.1.2 WATER 



For positive electrodes only the oxygen of the water molecule binds to the 
electrode metal atoms, since the strong electric field E will orient the large 
dipole fiw of the water molecule. The preferred position at positive potential 
bias is that of an inverted tripod, with two of the hydrogens pointing upwards. 
The adsorption fugacity of the 'oriented' water can be written as 

zw = Pw>^w9f (9) 
wherepvy = 3.345A~^ is the density of water and the adsorption parameter 

= e^^-^^(^'^-) AV^i^ = {i,- i,^-) (10) 

Here is the orientation parameter of the water molecules in contact 
with the electrode, pvy is the bulk density of water. We use the Gouy- 
Chapman theory to get the perpendicular component of the electric field at 
contact with the electrode 

Ik 

E^iA^pw) = (-)-^sinh(A^/'vK/2) where A^l^w = 38 MA(j){Volts) 

(11) 

A simple estimate of the orientation parameter of the water molecules near 
a charged electrode,(7^ can be obtained using a suitably adapted form of the 
mean spherical model [p2|, pSl El]: 



w _ V^f^jdwE^iAipw] 



A(2 + A) 



[1 - 1/8(A + l)r(T] (12) 



To conform to the high coupling- low concentration limit we use the EXP 
approximation 

^'^ " A(2 + A) exp{-(l/4)(A+ l)ra} (13) 

The parameters in this equation are k, the Debye Hueckel screening length, F 
is the MSA screening length, fiw, the dipole moment of water and a, the 
diameter of the water molecule (2.8 angstrom). The polarization parameter 
A is the MSA Wertheim parameter, obtained from 

4v^=A(l + A)2 (14) 
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ORIENTATIONAL ORDER 




Figure 2: Orientational order parameter. 



where ew is the dielectric constant of water, (78.4 at 298 "K), for which we 
get 

A = 2.652 



From these we get for the dipole model 
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eA(2 + A) 



sinh exp {-(1/4)(A + l)Ta} 



(15) 



or 
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0.0943 sinh (19.46A0vk) exp {-2A635Ta} (16) 

which is used in conjunction with Eq.(|l]) to compute the average orientation 
opf the water molecules. The results of this calculation are shown in Fig. 
2(1^), and the contribution to the differential capacitance (in region 3) are 
shown in Fig.3(|2]L2]). 



This procedure will be extended to include the full Yukagua model |25 



in the future. We do not expect however large differences since the main 
term in the re-orientation hamiltonian is the dipolar one, and the predicted 
capacitance for region 3, obtained from this model agrees reasonably well 
with experiment. The chain structure is shown in Fig.4 (|2rL^ ). 
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Figure 3: Orientational contribution to the capacitance / 
voltammogram. 



3 REGION 3: THE WATER CHAIN RE- 
GION 

The STM structures seen in these systems in the positive end of the voltam- 
mogram by various authors p^, |2^, |23], consists of elongated -\/3 x 



hexagons rather than regular hexagons. The elongated structures correspond 
to water coadsorption as was shown by Ito. More recent high resolution 
STM experiments P, ^ have been interpreted as being the result of hydro- 
gen bonded water chains. However, the structures that are proposed are 
unstable in the presence of the high field gradients |29| of the electric double 



layer, since the dipoles are alternated in order to keep the hydrogen bond 
straight. A different model is proposed in this paper in which all dipoles 
are polarized the same way. However this requires bending of the hydrogen 
bond. Our point here is that this is actually the case, even in normal ice 
Ih g, and very clearly in liquid water This model also explains why 



when the polarity is changed, those chains no longer exist and the regular 
a/S X hexagons are observed 0, ||]. The details of this model, as well as 



its consequences will be discussed in a separate paper [^, |T1] 



The formation of water hydrogen bonded chains occurs because the prop- 
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erly oriented water molecules form hydrogen bonds with the lone pair and 
one of the two hydrogens of the 'upper' plane. In principle any of the two 
can form the next bond, and therefore the 'water chain' chain will consist 
of units forming an angle of 142° to 152° instead of the straight 180°, and 
the 'free' hydrogen will be pointing to either side of the chain. If consecutive 
free hydrogens point to opposite sides of the chain, then the space between 
next nearest neighbors can be occupied by a (bi) sulfate ion. If the nearest 
neighbors free hydrogens point to the same side, then the chain will bend, 
and there will be a steric hindrance for the (bi) sulfate ion in the neighboring 
site. The model can be mapped into an equivalent one dimensional polymer 
with internal degrees of freedom. As such no first order transition occurs in 
the formation of these chains of elongatged hexagons. A detailed discussion 
of this model is left for a future publication 



3.1 Contribution to the voltammogram 



The main contributor to the current is the capacitive current due to the 
flipping of the water molecules, which is computed from the orientational 
fugacity Eqs.(|l],|16|) using a Langmiur type equation 



e 



W — T 



(17) 

where a is a constant that determines the position of the maximum of the 
voltammogram in region 3, and has been adjusted for the time being. The 
change in the orientation of the water molecules will produce a change in the 
capacitive current. 



jc = (iM) 



dtp ,,,,,, ^dO dih 
-r oc (M A)eC—-r- 
dt ^ ' ' ^dil) dt 



(18) 



where C is the integral capacitance, M/ A is the number of adsorption sites 
per area, e is the elementary charge, and C, is the partial charge of the adsor- 
bate. 



4 REGION 2: THE CP BUTTERFLY 

The classic shape of the cyclic voltammogram of the hydrogen UPD on a 
platinum single crystal with (111) orientation has been likened to the silhou- 
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Figure 5: The honeycomb bisulfate-water-hydronium lattice. It 
is commensurate with the adsorption sites of the Pt(lll) face. 
The © corresponds to a hydronium, the W to a water with the 
H pointing downwards, and the points are either an empty site 
or a bisulfate ion. 



ette of a butterfly, hence this particular CV profile is known as Clavilier's 
papillon{CP)^. This system has been studied experimentally by a large 
number of authors B, I 



34, 26, 27, 35, M, 37, 



The main point of our theory is the formation of an hexagonal two di- 
mensional honeycomb two dimensional 'ice' structure ( the HER( Hydrogen 
Evolution Intermediate) structure shown in Fig.5(^) that remains in place 
after the desorption of the bisulfate ion. The spike of the Clavilier papillon 
(CP, region 2)corresponds to the formation of an ordered, x y/S, com- 
mensurate phase 



-SOaH' 



-HoO + -mo' 



Lattice. 



(19) 



This is shown in Fig. 5 The 'wing' of the CP [42, H3, Mis associated with the 



desorption of (bi)sulfate, and corresponds also in our theory to the desorption 
of the bisulfate. However in our case the HER intermediate, the water- 
hydronium honeycomb (WHH) lattice stays intact, and is well represented 
by a Langmuir adsorption isotherm. This desorption mechanism is different 
from the hard hexagon order-disorder transition, recently proposed by Koper 
|5[, since in their case the desorption produces necessarily a disordered 



et al. 
lattice 
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The CO displacement experiments of Feliu and co-workers PBI, show 



very clearly that the 1/3 corresponding to the 'wing' of the butterfly is capac- 
itive and is due to the desorption of the bisulfate, which is what we propose 
in our theory. 



The detailed PCVM treatment will be published elsewhere [^. Some high- 
lights are: 

• Input parameters are the fugacity of 'up' water (W), the fugacity of 
'up' hydronium (H) and of the sulfate (S). However the formation of 
the honeycomb structure requires an equimolecular mixture of W and 
H, so that for the cluster theory we only need an effective adsorbate, 
which is one half of the W-H compound, the semi-hydronium (SH). In 
our model the assymptotic occupancy is 1/3 of W, 1/3 of H, or 2/3 for 
the SH complex, but the effective charge of the SH complex is only 1/2 

• the relevant parameters are the S-S interactions, the SH-SH=W-H in- 
teractions and the SH-S interactions. 



The result of the CVM are a series of sharp and soft transitions[^. In our 



case there is only one sharp transition, which corresponds the co-adsorption 
of bisulfate and hydronium and presumably, also the reorientation of water. 
The adsorption isotherms are represented by smoothened step functions: 
For the bisulfate 

es = hi + Erf[As{^ - ^s)]} + ^{-1 + Erf[-AsH{^ - ^sh)]} (20) 

D 15 

and for the semihyronium complex 

OsH = ^{1 + Erfl-AsHii^ - i^sn))]} (21) 

The positions and widths of the transitions are treated as adjustable pa- 
rameters, which however correspond the well defined physical parameters: 
^SH = 40 is the semihydronium transition width; 

'ipSH = 0.5041 is the position of the semihydronium transition, which occurs 
when there is simulataneous co-adsorption of bisulfate and semihydronium. 
As = 12 is the sulfate; 

ips = 0.44 is the position of the sulfate transition, which is determined by 
the bisulfate-bisulfate repulsive interaction. 

The bisulfate goes from a coverage of 1/5 to 1/3 , and the semihydronium 
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Hydronium-Bisulfate Coadsorption: Coverages 
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Figure 6: Semihydronium-bisulfate coadsorption transition 
erages. 
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Hydronium-Bisulfate Coadsorption: Current 
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Figure 7: Semihydronium-bisulfate coadsorption transition cur- 
rents. 
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from to 2/3. The exact values for the bisulfate coverages will depend on 
the other interaction parameters, and will be discussed in a more detailed 
PCVM analysis. The coverages for the semihydronium are fixed by the ex- 
perimental hydrogen yield. The results of the treatment for the adsorption 
isotherms are shown in Figs. 6(H) and 7(0). The agreement with the recent 
work of Kolics and Wieckowski|S] is qualitative: First the systems are not 
exactly the same, and in our case we predict a small glitch when the sul- 
fate is re-adsorbed. We remark also that it is very likely that bisulfate and 
perchlorate co-adsorb [| , since they are isostructural and isoelectronic. This 
point will be discussed in the future. 



THE UPD REGION 1: THE MODEL FOR 
THE HER INTERMEDIATE 



Region 1, which is called the H-UPD region psfhas been explored experimen- 



tally and several proposals have been made for the HER intermediate 



based on interpretations of the IR spectra. Our model is based on the 
stochiometry and the analysis of the phases of the voltammogram, the STM 
experiments and the recent radiotracer adsorption experiments, and is con- 
sistent with the experiment described by Peremans et al |^5[ Our HER model 



is a hydrogen bonded network of hexagonal rings, wwhich is a form of two 
dimensional ice. The chemical formula for a species is 

[{H,02)% (22) 

which forms a network of hexagonal rings. The explicit reaction taking place 
at the (111) face of the electrode is 

[{H^02y\i + 6e- = m2 + + (23) 

where the 6 hydrogen ions that are neutralized are those 'trapped' in the hol- 
low sites of the Pt(lll) surface. This means that the spectra observed in the 
SFG experiments of Tadjeddine are likey to correspond to Pt3 — H — — H-^/2 

^ Pechlorate is tetrahedral, the bond length is 1.325 A in the crystal. Bisulfate cannot 
be measured experimentally. The in vacuum structure has been calculated using state of 
the art DFT by S. Gaurei, is a slightly distorted tetrahedron with a S-0 distance of 1.43 
A. The major difference appears to be the dipole moment, which is zero for perchlorate 
and 3.80 D for bisulfate 
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(semi-hydronium) complexes. The 'upd' hydrogen is really a 'semihydro- 
nium' ion, with the hydrogen stuck in the hollow site of the platinum lattice, 
and the dipole of the water molecule properly oriented. Recent quantum 
mechanical DFT calculations support this model |^ . 



We assume that a 'surface state' is formed by the HER inermediate that 
discharges linearly with potential. Our model explains the 2/3 yield for the 
(111) face in a very natural way. The contribution of this region is only 
faradaic. 



jF = {M/A)e 



(24) 



where M is the number of adsorption sites per area A, e is the elementary 
charge, u is the electrovalence of the adsorbate in the bulk, ( is the par- 
tial charge of the adsorbate at the surface, and ip the potential, (z/ — () is 
the charge per adsorbate transferred to the surface. Fig. 8(|) shows the 
contributions to the faradaic current due to the processes of adsorption and 
discharge of the hydrogen adsorbate. Fig. 9(^ shows the rather nice agree- 
ment with the experiment ( courtesy of Dr. V. Climent) for a O.IM sulfuric 
acid solution. 



6 DISCUSSION OF RESULTS 

The main points in our work are that the the tetrahedral geometry of the 
water potential and the large dipole create two different water 'phases' for 
positively and negatively charged electrodes: Long chains at positive poten- 
tials and a honeycomb lattice for negative potentials. The formation of the 
honeycomb two dimensional ice phase is promoted by the bisulfate, or by 
other anions present, and the spike observed experimentally corresponds to 
this coadsorption. The honeycomb structure remains in place after the bisul- 
fate has left, and is responsible for the 2/3 hydrogen yield of the Pt(lll) face 
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Contributions to the total current 
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Figure 8: Contributions to the current. 
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Figure 9: Comparison between our model and the experimen- 
tal voltammogram.The experimental data is courtesy of Dr. V. 
Climent 



19 



G. Jerkiewicz, V. Climent,J.M. Orts and A. Wieckowski for the data files of 
voltammograms, and Prof. R. Car and Dr. N. Marzari for useful discussions. 



References 

[1] J. Clavilier, R. Faure, G. Guinet, and R. Durand, J. Electroanal. Chem., 
107, 205 (1980). 

[2] J. Clavilier, J. Electroanal. Chem., 107, 211 (1980). 

[3] V. Climent, personal communication. 

[4] D. A. Huckaby and L. Blum, J. Electroanal. Chem., 315 255 (1991). 

[5] D. A. Huckaby and L. Blum, in Microscopic Models of Elec- 
trode/Electrolyte Interfaces J.W. Halley and L. Blum, editors. The Elec- 
trochemical Society, Proceedings 93-5 232 (1993). 

[6] L. Blum and D. A. Huckaby, J. Electroanal. Chem. 375 69 (1994). 

[7] L.J. Wan, S. L. Yau and K. Itaya J. Phys. Chem.99 9507(1995). 

[8] L.J. Wan,T. Suzuki, K. Sashikata, J. Okada, J. Inukai and K. Itaya J. 
Electroanal. Chem.461 26 (2000). 

[9] J. G. Kim, J. B. Soriaga, G. Vigh and M. P. Soriaga, J. CoU. Interf, Sci. 
15 26 (2000). 

[10] J.M. Orts, L. Bum, D. A. Huckaby J.M. Feliu and A. Aldaz, Elec- 
trochim. Acta (in preparation). 

[11] L. Blum and J. M. Orts, in preparation, 

[12] M.L. Rosinberg, J.L. Lebowitz, and L. Blum, J. Stat. Phys., 44, 153 
(1986). 

[13] D.A. Huckaby and L. Blum, J. Chem. Phys., 92, 2646 (1990). 

[14] L. Blum and D.A. Huckaby J. Chem. Phys., 94, 6887 (1991). 

[15] E.A. Guggenheim and M.C. McGlashan, Proc. R. Soc. London, 206, 
335 (1951). 

20 



[16] G.M. Bell, J. Phys. C, 5, 889 (1972). 

[17] M. Shinmi and D.A. Huckaby, J. Chem. Phys., 84, 951 (1986). 

[18] D.A. Huckaby, M. Shinmi, M. Ausloos, and P. Clippe, J. Chem. Phys., 
84, 5090 (1986). 

[19] M. Shinmi and D.A. Huckaby, J. Phys. Soc. Jpn., 58, 571 (1989). 

[20] D.A. Huckaby, M.D. Legault, and L. Blum, J. Chem. Phys., 109, 3600 
(1998). 

[21] L. Blum, D.A. Huckaby, and M.D. Legault, J. Chem. Phys., , submitted 
(2001). 

[22] L. Blum and D. Henderson, J. Chem. Phys. , 74 1902 (1981). 

[23] L. Blum, F. Vericat and W.R. Fawcett, J. Chem. Phys., 96 3039 (1992). 

[24] L. Blum and W. R. Fawcett, J. Phys. Chem., 97 7185 (1993). 

[25] L. Blum, F. Vericat and L. Degreve , Physica A 265,396 (1999) 

[26] H.Ogasawara, Y. Sawatari, J. Inukai and M. Ito, J. Electroanal. Chem. 
358, 337 (1993) 

[27] Y. Sawatari, T. Sueoka, Y. Shingaya and M. Ito, Spectrochim. Acta, 
50A 1555 (1994) 

[28] A.M. Futnikov, U.Linke, U. Stimming and R. Vogel, Surf. Sci. 324, L343 
(1995). 

[29] N.F. Mott and J. Watts- Tobin, Electrochim. Acta 4 79 (1961). 
[30] R. Chidambaram, Acta Cryst. 14 467 (1961). 

[31] Y. Shingaya and M. Ito Interfacial Electrochemistry: Theory, Exper- 
iment and Applications A. Wieckowsky , Editor, M. Dekker, (1999), 
Chapter 17, p. 287. 

[32] A.S. Dakkouri and D.M. Kolb Interfacial Electrochemistry: Theory, Ex- 
periment and Applications A. Wieckowsky , Editor, M. Dekker, (1999), 
Chapter 10, p. 151. 

21 



[33] S. Motoo and N. Furuya, J. Electroanal. Chem. 172, 339 (1984). 

[34] F.T. Wagner and P.N. Ross Jr., J. Electroanal. Chem., 250, 301 (1988). 

[35] A.Peremans and A. Tadjeddine, Phys. Revs. Letters 73, 3010 (1994) 

[36] C. Nishihara and H. Nozoye, J. Electroanal. Chem. 379, 527 (1994). 

[37] R. Gomez, J.M. Feliu and A. Aldaz, Electrochimica Acta, 42, 1675 
(1997) 

[38] R. Comcz, A. Fernandez- Vega, J.M. Feliu and A. Aldaz, J. Phys. Chem., 
97, 4769 (1993) 

[39] J. Clavilier, R. Albalat, R. Gomez, J.M. Orts, J.M. Feliu and A. Aldaz,J. 
Electroanal. Chem. 330, 489 (1992). 

[40] A. Zolfaghari and G. Jerkiewicz, J. Electroanal. Chem., 420, 11 (1997). 

[41] A. Zolfaghari and G. Jerkiewicz, J. Electroanal. Chem., 422, 1 (1997). 

[42] K. Kunimatsu, M.G. Samant, and H. Seki, J. Electroanal. Chem., 272, 
185 (1989). 

[43] P.W. Faguy, N. Markovic, R.R. Adzic, C.A. Fierro, and E.B. Yeager, J. 
Electroanal. Chem., 289, 245 (1990). 

[44] A. Kohcs and A. Wieckowsky, J. Phys. Chem., , (2001) 

[45] M.T.M Koper and J.J. Lukkien, J. Electroanal. Chem.485 161 (2000) 

[46] J.M. Feliu, J.M. Orts, R. Gomez, A. Aldaz and J. Claviher, J. Elec- 
troanal. Chem.372 265 (1994) 

[47] B.Alvarez, V. Climent, J.M. Feliu, and A. Aldaz. Electrochem. Commu- 
nications2 427 (2000) 

[48] A. Bewick and A.M. Tuxford, Electroanal. Interf. Electrochem. 47 255 
(1973). 

[49] A. Bewick, K. Kunimatsu, J. Robinson and J.W. Russell, J. Electroanal. 
Chem.119 175 (1981). 



22 



[50] A. Bewick and J.W. Russell, J. Electroanal. Cliem.132 329 (1982). 

[51] N.Marzari and R.Car, personal communication. 



23 



